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ABSTRACT 

We study the mass distribution in six nearby [z < 0.06) relaxed Abell clusters of 
O ■ galaxies A0262, A0496, A1060, A2199, A3158 and A3558. Given the dominance of 

I dark matter in galaxy clusters we approximate their total density distribution by the 

NFW formula characterized by virial mass and concentration. We also assume that the 
' anisotropy of galactic orbits is reasonably well described by a constant and that galaxy 

^ \ distribution traces that of the total density. Using the velocity and position data for 

120-420 galaxies per cluster we calculate, after removal of interlopers, the profiles of 
the lowest-order even velocity moments, dispersion and kurtosis. We then reproduce 
the velocity moments by jointly fitting the moments to the solutions of the Jeans 
equations. Including the kurtosis in the analysis allows us to break the degeneracy 
I between the mass distribution and anisotropy and constrain the anisotropy as well 

f*^ , as the virial mass and concentration. The method is tested in detail on mock data 

extracted from A^-body simulations of dark matter haloes. We find that the best- 
fitting galactic orbits are remarkably close to isotropic in most clusters. Using the 
fitted pairs of mass and concentration parameters for the six clusters we conclude that 
O ' the trend of decreasing concentration for higher masses found in cosmological iV-body 

simulations is consistent with the data. By scaling the individual cluster data by mass 
C/3 . we combine them to create a composite cluster with 1465 galaxies and perform a 



m 
(N 



similar analysis on such sample. The estimated concentration parameter then lies in 
the range 1.5 < c < 14 and the anisotropy parameter in the range —1.1 < (3 < 0.5 at 
the 95 percent confidence level. 



' Key virords: galaxies: clusters: general - galaxies: clusters: individual: A0262, A0496, 

5^ , A1060, A2199, A3158, A3558 - galaxies: kinematics and dynamics - cosmology: dark 

matter 



1 INTRODUCTION 

Studies of galaxy kinematics in clusters remain a major tool 
in determining the mass distribution in these objects, com- 
plemented by methods based on the analysis of the hot X- 
ray gas and gravitational lensing. Due to a limited number 
of measured galaxy redshifts per cluster such analyses have 
been usually performed on composite clusters by combining 
data from many objects (e.g. Carlberg et al. 1997; van der 
Marel et al. 2000; Biviano & Girardi 2003; Mahdavi & GeUer 
2004; Katgert, Biviano & Mazure 2004; Biviano & Katgert 
2004; Goto 2005; Biviano & Salucci 2005). The normaliza- 
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tions needed in stacking the clusters together reduce however 
the number of parameters that can be estimated from the 
analysis. Besides, such studies are usually restricted to the 
analysis of velocity dispersion profile with the simplifying 
assumption of isotropic galactic orbits. 

In this study we attempt a kinematical analysis of six in- 
dividual clusters which we supplement in the end by a simi- 
lar procedure performed on a composite cluster created from 
the galaxies belonging to the six clusters. Our method re- 
lies on an extension of the usual Jeans formalism beyond the 
lowest-order velocity moment and including also the kurtosis 
of the line-of-sight velocity distribution (Lokas 2002; Lokas 
& Mamon 2003). The formalism has been successfully ap- 
plied to study the dark matter distribution in the Coma clus- 
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Table 1. Observational parameters of the clusters. 



assigned 


cluster 


RA 


Dec 


redshift 


number 


velocity dispersion 


kurtosis 


number 


name 


(J2000) 


{J2000) 


z 


of galaxies 


o-ios [km s~^] 


/^los 


1 


A0262 


0l'^52™50.4= 


+36°08'46" 


0.0163 


120 


527 ± 34 


2.63 ± 0.44 


2 


A0496 


04'^33™37.1= 


-13° 14' 46" 


0.0329 


270 


719 ± 31 


3.36 ± 0.30 


3 


A1060 


10'^36™51.3= 


-27°31'35" 


0.0126 


330 


696 ± 27 


2.68 ± 0.27 


4 


A2199 


16'^28™37.0" 


+39°31'28" 


0.0302 


180 


795 ± 42 


2.49 ± 0.36 


5 


A3158 


03'^42™39.6" 


-53°37'50" 


0.0597 


145 


970 ± 57 


2.58 ± 0.40 


6 


A3558 


13'^27™54.8= 


-31°29'32" 


0.0480 


420 


948 ± 33 


2.70 ± 0.24 



ter of galaxies by Lokas & Mamon (2003). It has been shown 
that, for a restricted class of dark matter distributions moti- 
vated by the results of cosmological A^'-body simulations, the 
joint analysis of velocity dispersion and kurtosis allows us to 
break the usual degeneracy between the mass distribution 
and velocity anisotropy and constrain the parameters of the 
dark matter profile. Recently we have tested the reliability of 
this approach against a series of A^'-body simulations (San- 
chis, Lokas & Mamon 2004) and also applied it to constrain 
the dark matter distribution in the Draco dwarf spheroidal 
galaxy (Lokas, Mamon & Prada 2005a). 

Here we further test the method on a different set of 
cosmological A''-body simulations by studying in detail the 
errors in the estimated parameters following from the sam- 
pling errors of velocity moments. We also introduce and test 
a new procedure of interloper removal which we then apply 
to the clusters. With a number of available galaxy redshifts 
per cluster much smaller than for Coma we had however to 
introduce a number of simplifications in our modelling com- 
pared to Lokas & Mamon (2003). We model the total mass 
distribution instead of only the dark matter component, we 
assume that galaxies trace the total mass distribution and 
we use the velocity data of all galaxies, not only ellipticals. 
On the other hand, we believe that the carefully selected 
clusters studied here are much more relaxed compared to 
Coma and therefore we avoid any uncertainties due to de- 
partures from dynamical equilibrium. 

The paper is organized as follows. In Section 2 we briefiy 
describe our data set. In Section 3 we summarize our method 
of data modelling including the removal of interlopers and 
fitting of velocity moments and test it on a sample of dark 
matter haloes extracted from a cosmological A'^-body simula- 
tion. The results for the six galaxy clusters and the compos- 
ite cluster are presented in Section 4. The discussion follows 
in Section 5. 



2 THE DATA 

We have searched the NASA/IPAC Extragalactic Database 
(NED) for nearby [z < 0.1) well-studied galaxy clusters with 
at least 120 galaxies within projected distance of about 2 
Mpc from the cluster centre and with cz velocities differing 
from the cluster mean by less than ±4000 km s~^. Among 
a few tens of clusters selected in this way we have chosen 
those which are likely to be relaxed judging by the regu- 
larity of their X-ray surface brightness maps obtained with 
Einstein (Jones & Forman 1999) and ASCA (Horner et al. 
2000) satellites. Next we have looked at the regularity of the 



diagrams showing the line-of-sight velocities of galaxies as a 
function of projected distance from the cluster centre. We 
have rejected those with irregular diagrams which may indi- 
cate merging or presence of neighbouring structures (Lokas 
et al. 2005b) for which our method of interloper removal 
does not work. 

Table lists the clusters we chose for the analysis to- 
gether with their positions and redshifts as given by NED. 
We also give the number of galaxies per cluster used for 
the calculation of the velocity moments (the numbers are 
for the final samples, after removal of interlopers). The last 
two columns give the values of the line-of-sight velocity mo- 
ments, dispersion and kurtosis calculated for all galaxies. 
The method of estimating these values is discussed in the 
next section. The list includes some very well known clus- 
ters. A0496 was studied in detail by Durret et al. (2000) 
who called it a prototype of a relaxed cluster. Indeed the 
cluster seems to have the most regular X-ray luminosity dis- 
tribution of those in our sample and a single central galaxy 
whose position coincides with the centre of the gas distribu- 
tion. The cluster A1060 (Hydra I, Fitchett & Merritt 1988), 
although quite similar to A0496 in many aspects, has two 
central galaxies and somewhat less uniform X-ray distribu- 
tion so it may have just reached equilibrium after a major 
merger. None of the clusters, however, is completely free of 
substructure (see e.g. Escalera et al. 1994) or neighbours 
(even A0496). A3558 (Shapley 8) is the richest member of 
the Shapley supercluster (Dantas et al. 1997; Bardelli et 
al. 1998) including a much smaller neighbour A3556 and 
a bigger, more distant cluster A3562. For this cluster we 
have restricted the analysis to distances where members of 
A3556 and A3562 are not likely to contaminate the sam- 
ple. A2199 seems to be quite a relaxed cluster but has a 
less massive neighbour A2197 and other groups (Rines et al. 
2001, 2002) which made us restrict the analysis to distances 
much smaller than the estimated virial radius. Also A0262, 
a member of the Perseus-Pisces supercluster, can be affected 
by neighbouring structures. 

It is generally believed that elliptical galaxies in clusters 
comprise a virialized, more relaxed subsample compared to 
(late) spirals which might be infalling into the cluster for the 
first time. It could therefore be desirable to restrict the anal- 
ysis to ellipticals as was done in the case of Coma cluster by 
Lokas & Mamon (2003) . However, for the present sample of 
clusters the morphological information is available only for 
a small fraction of galaxies making such an analysis impossi- 
ble. Given a larger number of galaxies per cluster one could 
also attempt to measure their surface density distribution 
reliably and use it as an input in the kinematical analysis. 
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With scarce samples presently available we are forced to as- 
sume that the number density of galaxies follows the total 
density distribution. 

For the kinematical analysis we have chosen as the cen- 
tres of the clusters their central cD galaxies which coincide 
with the centre of the X-ray surface brightness distribu- 
tion. In the case of the presence of two central galaxies (in 
A1060 and A3158) we have chosen as a centre the posi- 
tion of the one which is closer to the centre of the X-ray 
surface brightness distribution. The galaxy velocities have 
been transformed to the reference frame of the cluster and 
in order to calculate the distances within clusters we have 
transformed the cluster velocities to the reference frame of 
the cosmic microwave background. The concordance cosmo- 
logical model (ACDM) with parameters Qm ~ 0.3, 57a = 0.7 
and h = 0.7 is assumed throughout the paper. 



3 THE METHOD 



(6) 



whose value is 3 for a Gaussian distribution. 

Having measured the line-of-sight velocities and pro- 
jected positions for a number of tracer particles or galax- 
ies one can estimate the profiles of velocity dispersion and 
kurtosis by taking n tracer particles per bin and using the 
following estimators of variance and kurtosis 
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3.1 Overview of the method 

In this section we summarize our method of determining 
the mass distribution from velocity moments, as developed 
in Lokas (2002) and Lokas & Mamon (2003), and test it 
against mock data obtained from A^'-body simulations. The 
method relies on fitting the solutions of the Jeans equations 
for the second and fourth velocity moments to the profiles 
of the moments determined from the data. For the velocity 
dispersion projected along the line of sight, (J\os, the Jeans 
formalism gives (Binney & Mamon 1982) 



O-los(-R) = 



I{R) 



Vr2 - B? 



dr 



(1) 



where v(r) and /(-R) are the 3D and the surface distribution 
of the tracer as a function of a true (r) and projected (R) 
distance from the object centre respectively. The parameter 



/3= 1 



(2) 



describes a relation between the angular ob and radial a,. 
velocity dispersions and characterizes the anisotropy of the 
tracer orbits. We will assume it here to be constant with 
radius and consider —oo < /3 < 1 which covers all interesting 
possibilities from radial orbits (/3 = 1) to isotropy (/? = 0) 
and circular orbits (/3 — oo). For constant /3 the radial 
velocity dispersion, (7r(r) in equation Q is 



va^(r) — r 



'2/3 



d<E> 



r V—— dr 

dr 



(3) 



where $ is the gravitational potential. 

For the fourth projected velocity moment we have 



g(r,R,(3) dr 



liR) Jr - i?2 

where g{r,R,l3) = 1 - 2/?i?V^^ + /^(l + /9)-R''/(2r^) and 
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fVriP — const) = 3r 



'2/3 



2/3 2. , 

r uarir)—— dr 
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In the following we will rescale the fourth moment to obtain 
the line-of-sight or projected kurtosis 



is the mean of velocities in a bin. Since in the case of study of 
galaxy kinematics in clusters the number of galaxies usually 
does not exceed a few hundred, and for our least sampled 
cluster in TableQ(A0262) is as low as 120, in order to have 
at least 4 data points in each of the velocity dispersion or 
kurtosis profiles we need to adopt a rather low number of n = 
30 objects per bin. Using the Monte Carlo method described 
in the Appendix of Lokas & Mamon (2003) one can then 
construct unbiased and Gaussian-distributed estimators of 
line-of-sight velocity dispersion s and kurtosis-like variable 
k 



1/2 



1/10 



(10) 



(11) 



where 5* and K are given by equations Q-©- The factor n— 
1 in equation llOII is the well known correction for bias when 
estimating the sample variance, valid independently of the 
underlying distribution. In IIH the factor 3/2.68 corrects 
for the bias in the kurtosis estimate, i.e. unbiased estimate 
of kurtosis is K' = 37^/2.68, while the rather complicated 
function of K' assures that the sampling distribution of k is 
approximately Gaussian. We have also checked that defined 
in this way the measured velocity dispersion and kurtosis in 
a given bin are very weakly correlated so the data points can 
be fitted as indepen dent. The standard errors in the case of 
s are s/y/2{n — 1) while in the case of k are approximately 
0.02 (for n ~ 30). In the following we assign these sampling 
errors to our mock and real data points. 

We will further assume that the total density distribu- 
tion in the studied objects (simulated haloes or galaxy clus- 
ters) and the distribution of the tracer are well approximated 
up to the virial radius r„ by the NFW formula (Navarro, 
Frenk & White 1997) 



(12) 



q{s) ^ Ac c^g(c) 
Qco 3s (1 -f cs)2 

where s = r/r„, Qcfi is the present critical density, Ac = 
101.9 is the characteristic density parameter, c is the con- 
centration parameter and g{c) = [In(H-c) — c/(l-|-c)]~^. We 
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define the virial mass and radius as tliose witli mean den- 
sity Ac = 101.9 times tlie critical density according to the 
spherical collapse model for the standard ACDM cosmology 
(see Lokas & Hoffman 2001). The surface distribution of the 
tracer I{FC) following from 3D density profile 1121 can be 
found, together with other properties of the NFW haloes, in 
Lokas & Mamon (2001). 

Having estimated the velocity dispersion and kurtosis 
profiles from the measured positions and velocities we can 
fit these data with the solutions estimating three 

free parameters: the virial mass A/„, concentration c and 
anisotropy /3. In the remaining part of this section we apply 
the method to a set of mock data in order to assess its vi- 
ability in reproducing the properties of real galaxy clusters. 
Our approach here is similar to that of Sanchis et al. (2004) , 
but in addition we address the problem of interlopers and 
we estimate the errors in the parameters due to sampling 
errors. 



3.2 Removal of interlopers 

For this work we used the results of a cosmological dark mat- 
ter simulation described by Wojtak et al. (2005) . The simu- 
lation was performed within a box of size 150 h'^ Mpc as- 
suming the concordance cosmological model (ACDM) with 
parameters = 0.3, f^A = 0.7, h = 0.7 and erg = 0.9. We 
focused on 10 massive and isolated haloes extracted from the 
final output of the simulation whose properties are listed in 
Table 1 of Wojtak et al. (2005). In order to emulate the ob- 
servations we place an observer at a distance of 100 Mpc 
from a centre of a given halo so that he will be able to see it 
receding with a velocity of around 7000 km s"'^. We assume 
that the line of sight is parallel to the x, y or z axis of the 
simulation box, with respect to which the haloes should be 
oriented randomly. The observer is located far enough from 
the halo so that the cone of observation can be approximated 
by a cylinder. We then project all particle velocities along 
the line of sight and the distances on the surface of the sky 
restricting the observations to the circle of projected radius 
R = rv on the sky. Next we reject all particles with veloc- 
ities differing from the mean velocity of the halo by more 
than ±4000 km s~^, as we did for the real clusters. This 
cut-off corresponds to at least 4(Tios (with aios calculated for 
all particles in the halo) so it is not very restrictive. From 
the obtained sample of particles we randomly draw 300 out 
of about 10* per halo. 

The velocities and positions of 300 particles obtained 
in this way for halo 1 (of mass Af„ — 7.5 x 10^* Mq and 
virial radius = 2.3 Mpc) observed along x, y and z axis 
(upper, middle and lower row respectively) are shown in 
Fig. Q In the plots shown in the left column of the Figure 
the particles that happen to lie inside the virial radius of 
the halo (as verified using 3D information) were coded with 
filled circles, while those outside the virial radius are shown 
with open symbols. In the right panel the coding is similar 
but the particles were divided into those which are bound 
(have velocities smaller that the escape velocity, v < fesc) 
or unbound to the halo. The axes along which the particles 
were observed are marked in the lower left corner of each 
panel together with the labels showing whether the distinc- 
tion between the particles was made with the criterion of 
Tv or Vesc ■ In the upper left corners of the plots we give the 
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Figure 1. Line-of-sight velocities of 300 particles sampled from 
halo 1 as a function of projected distance from the centre of the 
halo. The upper, middle and lower row show observations per- 
formed along X, y and z axis of the simulation box respectively. 
Filled (empty) circles mark particles which actually reside inside 
(outside) r„ (left column) or which are bound (unbound) to the 
halo (right column) . The numbers in the upper left corner of each 
panel give the numbers of filled+empty circles. Solid lines show 
fitted ±3(Tios(i?) profiles separating the particles included in the 
analysis from those rejected as interlopers. 



numbers of particles fulfilling or not the criterion which sum 
to the total number of 300. 

Averaging over 10 haloes and 3 directions of observa- 
tion we find that of the total number of particles (300), 76 
percent reside inside rv, 87 percent inside 2 rv, while 92 per- 
cent are actually bound to its halo. In addition, the unbound 
particles are always a subsample of those with r > 2rv, i.e. 
there are no unbound particles inside 2r„. The particles that 
do not fulfill any of the criteria are more common at larger 
projected distances from the centre of the halo. The particles 
from outside r„ are obviously candidates for interlopers since 
they would not be used to estimate the density profile from 
the 3D information. However, as the cited numbers show, 
about half of them are actually close to the halo (within 
2 Tv) and most of them are bound to the halo and there- 
fore probably reasonably good tracers of the potential. This 
agrees with the recent studies based on A'^-body simulations 
(e.g. Klypin et al. 2003; Prada et al. 2005; Betancort-Rijo 
et al. 2005) which demonstrate that the virialized region ex- 
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tends somewhat beyond the virial radius, as we define it. 
Besides, the simulated dark matter haloes are usually not 
spherical and by imposing spherical symmetry in our defi- 
nition of virialized region we may in fact cut out particles 
which are actually members of the halo. Anyway, most of the 
candidate interlopers reside close to the mean velocity of the 
halo and only a few of them are true outliers which could 
significantly alter our estimates of the velocity moments. 

We proceed to remove these outliers in the following 
way. First we calculate the velocity dispersion profiles by 
binning the data (with 30 particles per bin) and assigning 
them sampling errors, as described in the previous subsec- 
tion. We then fit the data with solutions (0 assuming /3 = 
and adjusting Mv and c. Although the /3 values of dark mat- 
ter haloes are mildly radial (with mean /3 « 0.3, see e.g. 
Wojtak et al. 2005; Mamon & Lokas 2005) the actual values 
of the parameters are not very important at this stage since 
our purpose now is only to reproduce the shape of the dis- 
persion profile and this can be done well with 2 parameters 
instead of 3. If the velocity dispersion increases strongly and 
the fit goes to values of c < 1 (the NFW profile would not 
make sense since the scale radius would be larger than the 
virial radius) we keep c = 1 and adjust M„ and /3 instead. 
We then reject all particles lying outside the mean velocity 
±3(Tios(i?) where a\os{R) is the velocity dispersion profile 
obtained with our best-fitting parameters. The procedure is 
repeated until no more particles are removed. In each itera- 
tion we also calculate new estimate of the mean line-of-sight 
velocity of the particles with respect to which the rejection 
is performed. 

The ±3(7ios(i?) profiles from the last iterations are 
shown as solid lines in Fig. As we can see, the proce- 
dure removes the most obvious interlopers which affect the 
velocity moments most strongly. The mean number of par- 
ticles removed in all 30 experiments is 14 which corresponds 
to 61 percent of the mean number of unbound particles per 
halo. In rare cases where the initial velocity dispersion pro- 
file from the data is strongly increasing the first fit may not 
remove any particles and it is necessary to repeat the fitting 
for a smaller number of data points. In such cases it may also 
happen that a member particle from the centre of the halo is 
removed. The effectiveness of the method will be compared 
to other methods of interloper removal used in the literature 
in a forthcoming paper (Wojtak et al., in preparation). 



3.3 Fitting of velocity moments 

We now assess the viability of our method by jointly fitting 
both velocity dispersion and kurtosis profiles calculated from 
our mock data for the simulated dark matter haloes. After 
removal of interlopers we are left with 8-9 data points for 
both dispersion and kurtosis which is of the same order as 
what we have for real galaxy clusters. The quality of the ob- 
tained fits and magnitude of errors in the parameters should 
therefore also be similar. 

We jointly fit the mock data for velocity dispersion and 
kurtosis with the predictions from the Jeans formulae 
adjusting the three parameters Mv, c and /3. The parameters 
which minimize ^-r© shown in Fig. |5| as filled dots for 10 
haloes observed along x axis of the simulation box (for the 
observations along the y and z axis the results are similar). 
The la errors in the parameters due to the sampling errors. 




123456789 10 
halo 




123456789 10 
halo 



/3 



123456789 10 
halo 

Figure 2. Fitted parameters of the simulated dark matter haloes 
observed along the x axis of the simulation box. The three panels 
show from top to bottom the virial mass M„ (in units of 10^'' Mq), 
concentration c and the anisotropy parameter /3. The filled dots 
show the parameters estimated from the joint fitting of velocity 
dispersion and kurtosis data with la error bars. Crosses mark 
values estimated from the full 3D information. 
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were found by exploring the 3-dimensional confidence region 
in the — P, 13 — c and Mv — c parameter planes and find- 
ing probability contours corresponding to la, 2a and 3a i.e. 

= - Xmin = 3.53, 8.02, 14.2. The parameters of the 
haloes obtained from the full 3D information are marked 
as crosses. The masses and concentrations were taken from 
Table 1 of Wojtak et al. (2005) and the f} parameters were 
calculated from the P profiles of the haloes by taking aver- 
ages in 10 radial bins inside the virial radius. 

The discrepancies between the fitted and true values of 
the parameters may be due to non-sphericity, presence of 
substructure, projected outliers and departures from equi- 
librium e.g. in the form of streaming motions. However, as 
we can see, the true values of the parameters are almost al- 
ways within the estimated la error bars of the fitted values 
(except for the mass estimates for haloes 5, 7 and 9 where 
the fitted values are somewhat lower). We conclude from 
this analysis that the sampling errors are the main source of 
error in this method. 

We also note that the method might possibly work even 
better for real galaxies in clusters than for our randomly se- 
lected particles. Although it would be worthwhile to test 
the method on simulated galaxies, at the present stage of 
the simulations this would be reduced to using subhaloes 
detected with standard halo-finding techniques. The distri- 
butions of subhaloes both in space and velocity are known 
to be biased with respect to those of dark matter particles 
(Diemand et al. 2004) and probably still suffer from over- 
merging problem. The density distribution of subhaloes is 
flat in the centre while both the particles in simulated haloes 
and galaxies in clusters have cuspy profiles. It would be very 
difficult to disentangle the effects mentioned above from the 
uncertainties due to the use of subhaloes. 

Besides, simulations including baryons suggest that 
clusters tend to be more spherical than pure dark matter 
haloes (Kazantzidis et al. 2004; Basilakos et al. 2005) which 
would reduce the projection effects due to non-sphericity of 
the systems. We also believe that cluster-cluster mergers can 
be more easily detected in real clusters than in a single fi- 
nal output of an A^'-body simulation of dark matter haloes 
by studies of the X-ray emitting gas and so our sample of 
clusters is probably more relaxed than the sample of dark 
matter haloes we studied. In the application of the method 
to real galaxy clusters we will therefore neglect other sources 
of errors and estimate the uncertainties in the parameters 
only from the sampling errors of the velocity moments. 



4 RESULTS 

4.1 Removal of interlopers 

We now proceed to apply the method of joint fitting of ve- 
locity moments to the sample of six clusters listed in Table0 
We start by plotting in Fig. |21 the line-of-sight velocities of 
galaxies with respect to the cluster mean velocity as a func- 
tion of projected distance from the cluster centre. In order 
to separate the galaxies which will be used in the calcu- 
lation of velocity moments from the supposed interlopers 
we apply the procedure for interloper removal described in 
the previous section, exactly as we did for simulated data 
with the additional assumption that the tracer density is 



3 - 
2 - 
1 



-1 

-2 
-3 - 
-4 



0.5 1 1.5 2 2.5 
R [Mpc] 




0.5 1 1.5 2 2.5 
R [Mpc] 



3 - 
2 




-1 

-2 - 
-3 - 
-4 



A3158 



J • _ ** ■ 



1 




0.5 1 1.5 


2 2.5 




R [Mpc] 














A2199 _ 


1. 






% 














* 













0.5 1 1.5 


2 2.5 


J 


R [Mpc] 





A3558 



/7 



0.5 1 1.5 2 2.5 
R [Mpc] 



0.5 1 1.5 2 2.5 
R [Mpc] 



Figure 3. Line-of-sight velocities of galaxies with respect to the 
cluster mean velocity as a function of projected distance from the 
cluster centre for the six clusters listed in Table 1. Solid lines show 
fitted ±3(Tios(R) profiles separating the galaxies included in the 
analysis (filled circles) from those rejected as interlopers (empty 
circles). 



proportional to the total mass density (we do not infer the 
tracer density from the observed surface number density of 
galaxies, because the latter would be very uncertain for such 
small samples and probably suffers from incompleteness that 
varies with projected radius). 

Since we do not a priori know the virial radius of the 
cluster, after estimating the virial mass in each iteration we 
check whether all fitted data points lie inside the virial ra- 
dius or an additional data point could be included and adjust 
the number of points accordingly. The ±3aioa{R) profiles ob- 
tained in the final iteration of the procedure are shown as 
solid lines in Fig.|3] As can be read from the Figure the num- 
ber of rejected galaxies is between (for A2199) and 4 (for 
A3558), much lower than the number of particles rejected 
during the application of the method to simulated haloes. 
The rejected galaxy in the centre of A3558 lies inside the fi- 
nal ±3(Tioa(i?) lines but was removed in the earlier iteration 
of the procedure. 

Having found the final sample of galaxies we divide the 
data into radial bins of 30 galaxies (except for A3158 which 
has 29 galaxies per bin) and calculate the velocity dispersion 
and kurtosis profiles which are shown in Fig. 2] and |^ with 
la sampling errors. We see that the profiles are similar to 
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Figure 4. Line-of-sight velocity dispersion as a function of pro- 
jected distance from the cluster centre for the six clusters. Solid 
lines show the best-fitting profiles with parameters listed in Ta- 
ble 
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Figure 5. Line-of-sight kurtosis variables as function of projected 
distance from the cluster centre for the six clusters. Solid lines 
show the best-fitting profiles with parameters listed in Table 1^ 



those characteristic for objects with NFW-like density pro- 
files and orbits close to isotropic: the global trend is that 
both profiles slightly decrease with radius (see Lokas & Ma- 
mon 2001; Sanchis et al. 2004). A3558 shows the most vari- 
able velocity dispersion profile with one discrepant point at 
0.7 Mpc (which is due to the single galaxy with discrepant 
velocity present in this bin - see Fig. I^J and a secondary 
increase at about 1.8 Mpc while its kurtosis profile remains 
rather uniform. On the other hand A0496, believed to be a 
very relaxed cluster, has the most variable kurtosis profile. 
In an attempt to verify whether this variability may be due 
to departures from equilibrium e.g. in the form of the pres- 
ence of infalling groups of galaxies we have also looked at 
the mean line-of-sight velocity profiles with respect to the 
cluster mean. They do not however depart strongly from 
zero, typically remaining in each bin within 0.2K, (the cir- 
cular velocity at the virial radius as determined in the next 
subsection). 

4.2 Estimated parameters of the clusters 

In the following analysis we assume that the galaxy distri- 
bution follows that of the total mass distribution. Although 
this assumption is far from obvious there is evidence from 
studies of clusters (e.g. Carlberg et al. 1997; Lokas & Ma- 



Table 2. Fitted parameters of the clusters. 
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mon 2003; Biviano & Girardi 2003) that both number den- 
sity and luminosity density of galaxies in clusters are cuspy 
and therefore can be quite well approximated by the NEW 
profile (contrary to the distribution of subhaloes in the dark 
matter simulations, which probably have a shallow core, see 
Diemand et al. 2004). The concentration of these distribu- 
tions does not have to be the same as that of the total mass, 
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Figure 7. The mass-concentration relation. The filled circles 
sliow the best-fitting M„-c pairs for the six clusters with Icr error 
bars. M„ is given in units of 10^^ Mq. The solid line plots the 
prediction from the A^-body simulations by Bullock et al (2001) 
while the dashed line shows the best fit to the data points. 



but with numbers of galaxies as low as 120 per cluster we are 
not able to reliably estimate the distribution of the tracer. 

We jointly fitted the data for velocity dispersion and 
kurtosis with the predictions from the Jeans formulae Q- 
^ by minimizing and adjusting the three parameters 
M„, c and /3. The best-fitting parameters for the six clus- 
ters are shown in Fig. |S| as filled dots together with la 
error bars. The exact values are listed in Table |5| where 
the last column shows also the goodness of fit measure, 
X^/N. (The virial radii corresponding to the virial masses 
are r-„ = 1.2[M„/(10^*Mq)]^''^ Mpc for our adopted cos- 
mological model.) The errors in the parameters due to the 
sampling errors were found by exploring the 3-dimensionaI 
confidence region in the M„ — P, j3 — c and A4v — c param- 
eter planes and finding probability contours corresponding 
to Icr, 2a and 3a i.e. Ax^ = - Xmin = 3.53, 8.02, 14.2. 
The best-fitting velocity dispersion and kurtosis profiles are 
shown in Fig. 01 and |S] as solid lines. 

Interestingly, A3558 has the lowest concentration 
among all studied clusters. Although the rather large er- 
ror bars prevent us from concluding too much, this is ex- 
pected since A3558 is in the Shapley supercluster and there- 
fore probably in the early stage of cluster-cluster mergers 
which destroy the inner cusp. It is also worth noting that 
the highest value of anisotropy {l3 — 0.26) is obtained for 
A0496 which has the highest value of line-of-sight kurtosis 
(see Table 1) as expected for non-rotating systems (Merritt 
1987). 



Figure 6. Fitted parameters of the clusters. The three panels 
show from top to bottom the virial mass M„ (in units of lO^^M©), 
concentration c and the anisotropy parameter /3 with Icr error 
bars. The numbering of clusters is the same as in Table HI and 1^ 



4.3 The mass-concentration relation 

The relation between the virial mass and concentration of 
dark matter haloes is a well established result of A'^-body 
simulations and has been studied by many authors (e.g. 
NFW; Bullock et al. 2001; Dolag et al. 2004). In this subsec- 
tion we address a question whether our best-fitting param- 
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eters of the clusters agree with the trend of concentration 
decreasing with mass as found in A''-body simulations. We 
will compare our results to those of Bullock et al. (2001) 
because they used the same definition of the virial radius. 

In Fig. |7| we plot estimated concentration c as a func- 
tion of cluster virial mass M„ (in units of IO^Mq) for our 
six clusters. The Icr errors following from sampling errors 
of velocity moments were assigned to the points. The solid 
line shows the approximation of the c(M„) relation calcu- 
lated from the toy model proposed by Bullock et al. (2001) 
which reproduces well the properties of a large sample of 
haloes found in their A^-body simulations. The predictions 
were obtained for the standard cosmological model as in our 
simulations and z = from their equations (9)-(13) with 
parameters F = 0.001 and K = 3.0 as advertised for masses 
M > lO^*/i-^M0. 

Although the errors in the estimated parameters are 
quite large it is clear from Fig. |71 that the data for clusters 
are consistent with the c(M) relation found in A^'-body sim- 
ulations. To make this statement more quantitative we have 
fitted to the data (neglecting the errors in mass) a linear re- 
lation of the form logc = a log[M„/(10"MQ)] -|-6 and found 
the best-fitting parameters a — — 0.6± 1.3 and b = 1.3 ± 1.1 
(at 68 percent confidence level) . This best-fitting relation is 
shown in Fig. |7| as a dashed line. Therefore, although the 
best-fitting slope is negative, the data are also consistent at 
Icr level with constant concentration or even concentration 
increasing with mass. 

4.4 The composite cluster 

The low number of galaxies per cluster in our sample results 
in rather large sampling errors in the measured velocity mo- 
ments and therefore also large errors in the estimated pa- 
rameters of the clusters. In order to reduce the uncertainties 
and study a typical cluster we have combined the position 
and velocity data for our six objects to create a composite 
cluster. We have normalized the galaxy distances by the esti- 
mated virial radius of their cluster, and the velocities by the 
circular velocity at the virial radius Vv = GMv/tv In this 
way we make the velocity moments independent of the virial 
mass. The normalized velocities and velocity moments (cal- 
culated with 121 galaxies per bin) for the composite cluster 
made of 1465 galaxies are shown in Fig. |H| The kurtosis-like 
variable is now k = [log (3i^/2.90)]^/^° where the coefficient 
2.90 was adjusted to the number of galaxies per bin. The 
combined samples of galaxies were those after removal of 
interlopers in each cluster, but we perform the procedure 
(fitting velocity dispersion profile and adjusting concentra- 
tion while keeping /3 = 0) again on the total sample and 
it turns out that two more interlopers have to be removed. 
The two galaxies have been marked as before with open cir- 
cles in the upper panel of Fig. |H| The two solid lines in this 
panel plot the last iteration of ±3cric,s(7?) profiles separating 
the interlopers from the galaxies included as members of the 
cluster. 

With the mass dependence factored out, the velocity 
moments of the composite cluster depend only on two pa- 
rameters: concentration c and anisotropy (5. We adjust the 
parameters first by fitting only the velocity dispersion pro- 
file. The best fit is then at c = 2.1 and 13 = 0.34 with 
Xmin/-^ ~ 6.4/9. The la, 2a and 3cr probability contours 
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Figure 8. Upper panel: line-of-sight velocities of 1465 galaxies of 
the composite cluster normalized to Vv as a function of projected 
distance in units of the corresponding virial radius. Solid lines 
show fitted ±3(Tios(iJ) profiles separating the galaxies included 
in the analysis (filled circles) from those rejected as interlopers 
(empty circles). Middle and lower panel: line-of-sight velocity dis- 
persion and kurtosis profiles for the composite cluster. Solid and 
dashed lines show the best-fitting profiles with parameters listed 
in the corners of the panels. 
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Figure 9. The Icr, 2cf and 3(t probability contours in tlie c — /3 
parameter plane obtained from fitting the velocity dispersion only 
(upper panel) and both velocity dispersion and kurtosis (lower 
panel) of the composite cluster. The dots mark the best-fitting 
parameters. 



corresponding to Ax^ = X"^ ~ Xmin = 2.30,6.17,11.8 are 
shown in the upper panel of Fig. El with the dot mark- 
ing the pair of best-fitting parameters. We can see that 
the Icr confidence region is not restricted to the vicinity 
of the best fit but appears also at more negative /? and 
larger c, towards the second local minimum at c = 14.1 
and P = -2.1 with xV^ = 8.1/9. Fitting both the ve- 
locity dispersion and kurtosis breaks this degeneracy but 
only to some extent: the best fit is now at c = 2.2 and 
/3 = 0.35 with Xmin/^ = 11-0/20, but another local mini- 
mum of x^/N = 12.9/20 is found at c = 10.1 and /3 = -0.35. 

The presence of the two minima can be traced to the 
parameters estimated for each cluster separately in the last 
Section. In our sample we had three clusters with concentra- 
tion around 3 and three clusters with concentrations of the 
order of 10. Those more concentrated also had lower mean 
value of j3. Indeed, performing similar fitting for composite 
clusters made separately of these two groups we find single 
minima at c = 1.8 and /3 = 0.2 for low-concentration clus- 



ters while c = 12.7 and P — —0.5 for high-concentration 
clusters. 

The profiles of the normalized moments for the best- 
fitting parameters are plotted as solid lines in the middle 
and lower panel of Fig.|H| The solutions corresponding to the 
second local x^ minimum are plotted as dashed lines. Recall, 
that these are best-fitting profiles obtained from the joint fit 
of dispersion and kurtosis and not to each of them separately 
- this is why the fitted kurtosis profile does not match the 
corresponding data perfectly. Actually the best fit to the 
kurtosis data alone would have (3 closer to zero (isotropic 
orbits), but the velocity dispersion forces (3 towards more 
radial or more tangential values. Therefore, to summarize 
our results for the composite cluster, we find 1.5 < c < 14 
and —1.1 < P < 0.5 at the 95 percent confidence level. 



5 DISCUSSION 

We performed a joint analysis of velocity dispersion and kur- 
tosis profiles of galaxies in six nearby relaxed galaxy clusters 
estimating two global parameters of their mass distribution 
and the anisotropy of galaxy orbits. The method of joint 
fitting the moments was tested on simulated dark matter 
haloes and a new procedure for interloper removal was ap- 
plied. With the presently available number of galaxy red- 
shifts per cluster the errors in the estimated parameters are 
still large. The estimates of mass have an error from 20 per- 
cent in the case of A0496 to 50 percent in the case of A3158 
and are in rough agreement with previous estimates (e.g. 
Girardi et al. 1998; Rines et al. 2003). The uncertainties in 
the estimates of concentration are even larger. The fitted pa- 
rameters are consistent with the mass-concentration relation 
found in A''-body simulations. 

The method has been first applied a few years ago to the 
Coma cluster (Lokas & Mamon 2003) for which the largest 
number of galaxy redshifts is available. This larger amount 
of data allowed for a more detailed analysis: the dark mat- 
ter distribution could be modelled separately from stars and 
gas, the luminosity density profile could be used as a tracer 
distribution without assuming that galaxies trace the total 
or dark matter distribution, galaxies identified as members 
of binaries could be removed from the sample and the anal- 
ysis could be restricted to early type galaxies which are be- 
lieved to be more relaxed. For the present sample of galaxy 
clusters we had smaller numbers of galaxies and had to sim- 
plify the analysis. However, we believe that the clusters dis- 
cussed here are really relaxed which appears not to be the 
case for Coma (which is probably the product of a recent 
merger as suggested by its perturbed gas distribution and 
the presence of two cD galaxies). In addition, in the present 
study we have improved the method by introducing a reli- 
able new procedure for the removal of interlopers. 

Inclusion of kurtosis in the analysis allowed us to con- 
strain the anisotropy and thereby break the degeneracy be- 
tween the mass distribution and anisotropy present in the 
analysis of velocity dispersion. This led us to expect that 
very tight constraints on the concentration parameter and 
anisotropy could be obtained if the analysis was performed 
on a composite cluster made of galaxies belonging to all six 
clusters. Although the anisotropy was indeed constrained 
also in this case some degeneracy still persists between c and 
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/3 and almost equally good fits can be found for very different 
pairs of the two parameters. The degeneracy may however 
also be due to the rather simple form of anisotropy which 
we assumed to be constant with radius. Although isotropic 
orbits {/3 = 0) can very well be a realistic case for virialized, 
mostly elliptical galaxies in clusters, our samples contained 
also spirals which probably are on more radial orbits (Bi- 
viano & Katgert 2004). While isotropy might be preserved 
in the centre, our composite cluster may be better described 
by an anisotropy parameter increasing with radius. This pos- 
sibility will be studied elsewhere. 

The critical factor in studies of galaxy kinematics in 
clusters is the number of available redshift measurements. 
Since the number of such measurements typically does not 
exceed a few hundred, while the errors in the estimated 
parameters are mainly due to sampling errors, as we have 
shown, the future of such studies will probably still be in an- 
alyzing composite clusters. This does not allow however to 
study relations between individual objects, like the mass- 
concentration or mass-temperature relation. Those have 
been reserved, till recently, to studies based on X-ray obser- 
vations. Although much tighter constraints e.g. on the con- 
centration parameter can be obtained in this case (Pointe- 
couteau, Arnaud & Pratt 2005), these studies rely on as- 
sumptions like the hydrostatic equilibrium and rather un- 
certain temperature profiles of the X-ray emitting gas. They 
should therefore be complemented by extensive analysis of 
galaxy kinematics in clusters. This may soon become possi- 
ble with systematic photometric and spectroscopic surveys 
like WINGS (Fasano et al. 2005). 
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